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ABSTRACT 

The problem of predicting the motion, of a space 
vehicle in the vicinity of & Trojan point in cis lunar 
space is considered. Specifically, it is desired to 
know if the vehicle will remain in this vicinity, or 
diverge from the Tirol an point. If the equations of motion 
of the vehicle are linearised at tha Trojan point, the solu- 
tion of these equations indicates that the vehicle will 
aaevaca a periodic ^notion about the Trojan point. Fuodsraan- 
tally, the problem reduces to determining the validity, or 
rang© of validity of those solution 3 » 

Various methods are discussed along with their aseo- 
elated problem*. Further investigation is proposed along 
the lines of identifying and reducing computational errors. 
Also it la proposed to investigate the significance of 
linearising the equations o£ motion at points displaced 
from the Trojan point* 
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OBJECT 



It is dsoired to study the orbital stability of a 
space vehicle at a Trojan point 1b cis lunar space, i.e. ; 
where the earthy moon and vehicle form an equilateral 
triangle in the orbital plane of the earth and moon. 2i 
© vehicle originally placed as: the Irojsn point is indead 
in a stable orbit relative to the Trojan point, it should 
then be determined what order of inaccuracies of injec- 
tion can fee tolerated to still establish a stable orbit. 

This assumes shat there is an operational require- 
mens for a space vehicle of this sort which maintains its 
relative position with respect to the earth and moon. The 
requirement sight be for a space navigation or ecs'tusica- 
tions satellite, or it might fee for a rep lenistsaenf: station 
for g manned deep space probe. However, the specific 
advantages of a Trojan point orbit war© not investigated 
further. 
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On© solution of fcfeo classical pj.obloi of fibres bodies 
is the equilateral itri&n gle solution* diseoversd fey Lcgzmgv 
in 1772 y, where the three bodie-- are each equidistant frosa 

the other tss© snd* in the sboenee of eternal perturbation^ 

* 

s»inta£m. this relative poc.it ion. In the restricted eh.. s-d~ 

body problem %ih«re the wet** of the third body io si-tall tuc 
not to influence tin sioticn of the ether two feodion, tKe fcwu 
equilateral t%i-tiigl3 points :os the pl-steia of rotation of the 
two taajor beriisv ere celled Trojan points • The Trofjjs 
asteroids,, the fir »t of which u-tes dioccveted in 19o6, he.s& 
been foond to retain in. the vicinity of the* quilateral 
fcrianjjia points with respect t<* Jupiter -i>«J the euu. 

The tezs* orbital stability, if- used in the se^-j© tUat 
a body in .a stable orbit in the vicinity of .-; Trojan point 
re,*felns indefinitely in this vi»vinity s i*e . 9 tiha&g Is -ri, 
•upper bound t<* its dist&ace fires t/he f &ojw po^rit -*ifcieb 
does mot iaeresifce with tises. 

If thie definition is .«ce;-pt*ed li to tolly w@ ssmst 
alse Ifeifc the aw^it’jdc ••*>£ vhr liy r fe-vr .d. mince t -iy 
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ci® lunar vehicle with less th*m escape velocity will regain 
in cialunar space, barring collisions with another body, end 
will haves ass upper bound to its displacement from any point 
in the earth-moon coordinate system. 4 convenient raetr lo- 
tion cn the rigorous definition or stability is not iassediately 
apparent. Hos'/ever, the distinction la perhaps overly meticu- 
lous. We are concerned with a vehicle which remain® in the 
vicinity of & Trojan point and, dug to initial displacement 
fre® the Trojan point, initial mon-sera relative velocity 
with respect to the Trojan point, or partubntions external 
to the system, appears to oscillate about the Trojan point. 

If it continues to do this, it is stable^ if it diverges free 
the Trojan point, it is unstable. 

Historically, the Trojan points have been considered sta- 
ble for "smll disturbances" based on a linearisation of the 
equations of motion of the third body at the Trojan point. 
However, the stability of a body ia eha actual non- linear system 
has not been conclusively determined. In fact, there are 
arguments indicating that, even granted ersst Im owledge of the 
physical constants and unlimited precision and speed of digital 
computation, the stability can still cot be established with 
mthematie&l rigor. 

The existence of the Trojan asteroids is not as good an 
argiEsant for stability as it seszno at first thought. According 
to the solution of their linearised equations of motion, the 
periods of their two modes of vibration are 11.86 years and 
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467 years. Tha corresponding periods in aiQhmm: space are 
28.6 days and 91.7 days. Thus the 54 years available for 
observation of the Trojan asteroids corresponds, in one 
sense, to only about XI days observation of a eislunar 
Trojan point satellite. 

We would like to know specifically the ranges of initial 
conditions and magnitudes of perturbations for which a Trojan 
point vehicle would be stable. However, the study described 
here does not reveal such a criterion. Rather, It seesss we 
fsusfc be content with specifying s else limit during which the 
vehicle will not exceed a specified displacement fre-a the Tro- 
jan point. The intent of this investigation was to provide 
answers to questions of the sort, ”Whafc are the guidance 
accuracy requirements for establishing a dolman Trojan point 
spas© station so remain within 20,00'.) miles of the Trojan point 
for at least lo years?” Again, it must be admitted that the 
answer is not obtained In the present investigation. However, 
the directions toward which such answers lie way be Indicated. 



CHAPTER 2 



RESTRXCTtC^S CM THE PROBLEM 

~rar »i nn « «r. — «>■>■ 

The problem is considered to he th$£ of the restricted 
three-body problem, where the mass of the third feody^ the 
vehicle, is negligible compared with the other two bodies. 
The vehicle trill then have s*o effect on the motion of the 
other two bodies. This is obviously valid for the problem 
of earth, moon and a vehicle launched from the earth. 

Another restriction is the elimination of all pertur- 
bative forces ©asternal to the earth-soon systems. This is 
equivalent to considering the berycentsr of the earth-moon 
syoteai as an inertial point, and neglecting the curvature 
and velocity variations of its orbital motion. Tha pertur- 
bative effect of the aun is perhaps essential in the solu- 
tion of the physical problem. However, it con later be 
added to the differential equations as a forcing function 
when the simpler system is more fully understood. Other 
perturbative forces such as that due to Jupiter, tidal 
forces, earth ©b lateness and radiation pressure aro probably 
Indeed negligible. 

Also the orbit of the Eason about the earth, or, racre 
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precisely, about tbs barycenter of: eh© earth-aiacm systesa, 
is considered to be circular. This effect can also be 
added later if desired. 

Errors due to these als?plifyiws seoussptions will be 
discussed further in the text . Their effect i© also con- 
sidered by Bucheisj (1) . 
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the htmA^rm zozms 



3.1 The Geooas tr y of a Space V e Me Igt jA tCrac t od by Two Bodie s 

In Fig. 3. l p the vehicle , S 2 C , is attracted fey the two 
celestial bodies, ®j and %>, which revolve about their i&eo 
tial fearycenter, B. In eloltmr space, ss-, and ea^ are the 
earth and moon., ©sd, for convenience, they will bo s® neaacc. 
The mass of tha vehicle is considered s*a!2 eaov^ia to have 
negligible effect on the motion of the earth and moon. 

To simplify the equations , all distances are aossalised 
with respect to the earth-soon distance. Mefcaneea faraaa 
earth and moon to their fearyceater are then - end ft respse® 
tively -dbszo 

e ^2 

<o«j* u3 c^ssscji:* 

A E5j 
6 -ir ft C3 1 



by 



The grasvifc&t i&att foreoo on the vehicle s 
'£-* gad T« which are vector specific fore sc. 



re represents ;* 
or force per 



unit K-aso of the vehicle. 








> 




(3.H 



> 



( 



* 







( 



c. 




Fig, 3.1= Geometry of & Space Vehicle Attached by Two Bodies 



(3.2) 
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It is cosrvo&ient to choose a rotating coordinate 
system in which the earth and moon appear stationary. 

The equations of asotlon in this coordinate system oust 

I 

then include centripetal and Coriolis accelerations 
resulting from the coordinate transformation. It is also 
convenient to consider the centripetal acceleration of the 
vehicle due to coordinate transformation as a centrifugal 
specific force, T , along the radius vector from ths bary~ 

■w 

center. 

£ c 0 rfi" 



where SI la the angular velocity of the coordinate system 
with respect to inertial space. Then free Kepler’s third 
lot? 

q2 ^ G (djd'3ig) 

♦Suy \^xam*xmcnas~ 

(l) i (3-3) 

£ n = G (ta^*S?ag) r 



3.2 The Existence of Libration Points 

The following development, due to Klessperer end 
Bsnedikt, (2) sham the insistence of five points in the 
plane of the moon’s orbit at which a vehicle with zero 
relative velocity in the rotating coordinate system 
© 2 £p©ri< 2 EC 0 O ser© resultant force. These points were 
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first discovered by L&graa. ge in 1772 as particular 
solutions of the three-body problem. Their existence 
could bo proved with less effort by simple substitu- 
tion of their coordinates into the equations of motion, 
but the following development also shows that these 
are tha only five points at which this condition occurs. 

Consider the vehicle to be at rest in the rotating 
system* With sero relative velocity* the Coriolis 
acceleration is aero and* if the vehicle is in equili- 
brium 



Consider first the circiaafereatisl equilibriisa of 
the vehicle* or the suKsaation of force components per- 
pcsndiculer Sc the radius vector from she b&rycenter. 



One locus of circumferential equilibrium is the straight 
line through end where 



£.j 4* Tg 4* T c « 0 




f « 0 (3*4) 




1 



2 




Also * by the law ©£ sines 



sa -a in iff ~ C P) » sis *0 



sin ^ 

*? 

Substituting for sin /&&&& sin f J is Eq. 3.4 
(€ aisa^ ) ^ g C?aia 0 ) 

■ n» ■ H H. . V SC 5 iw « rrt»l ■« Vx >» f«rn » 

-3 

r l r 2 

*«d tbd other locus of circnjafereiafcial equilibrium Is 
the line 

r l s r 2 




Hesit consider fcfea radial equilibrium of the vehicle, 
or the ousssatioa of force components parallel to the radius 
vector free the fearyeest&r . 



CQ3 JC » fg COS/^4* f c a o 



Itt-idtS*} \ 
is—feL 
2 



COS 



r 






2 
1 2 

cos 2f ~ 



co© ^4- rG (mj^sig) «* 0 



cos 4* 4* r w C 

■nap— 



Finding cos X and eos^by the Xss-? of cosines 
cs 4* r| «• 2rr^ cos 



(3.5) 



cos^® r -1* r. 



1 ~^r 



2rr-, 



.2 



2 



«» ff r: 



.2 



2 




lx 





V 



















COB $ w 



"2 



Zrr, 



Substituting for c os Z. ©d cos ^ in Sq. 3. 



4- (jf- + r i - * 2 1 - e 

r[ V2r l 1r ^ 



«2 



~ **’ r 2 Y& ~ 1 + * 

2 



.<*- I 2^ ' "'2rr 



0 



3 3 

Multiplying by »2&r£ ~2 to fcb® fractions 

^?r 2 r| -Mr^ Sr2r i * er i r 2 “ e ^ 2?r i " 2r2r l r 2 3 “ J 

<3. 6) 

To eliminate r £ross the equation, y is found is* teres 

of r-^ and by the las* of cosines. 

■% 

r 2 » € 2 + r? “ 2er„ coa 

and also 

r| « 1 -3* r 2 - 2r^ cos 



(3.7) 

(3.8) 



Substituting for^ cos^f from Eq. 3*7 into Eq. 3.8 
r 2 » e 2 * r 2 - s (Xi*r|-r 2 ) « e 2 -e * r| (1-e) 4* e® 2 

r 2 » «(i-e) + y | (I~e) "3* «r 2 » 4- e^ 2 

2 

Substituting for r in Eq. 3.6 and cabining terns ^ the locus 
for radial equilibrium is obtained 



2 J 3 
1T~ 



5 3 , 2 3 5 9 3 3 5 .J 1+s 32 1+6 _2 3 

r l r 2 + 7T r i r 2 * 2l X r 2 ■ r l* f 2' "S - r l*2 " -r- r l r 2 ’• 



3 3 

+r^ 4- r 2 « 0 



The equilibrium loci are shown in Fig. 3.2. The five 
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intersections of the two sets of curves are eha only points 
at which a body would rssnin at equilibrium. Those five 
points ©re statically stable in the rotating coordinate 
system, but the development yields ao dynamic information. 
To specify the behavior of a vehicle ia the vicinity of a 
Vibration point with deviations in position and non- aero 
relative velocity, the equations of motion of the vehicle 
must be investigated. 




Fig. 



3.2c 



LiOCfjcioti of the bibration Foists 
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CHAPTER 4 



THE EQUATICHS OF MOTION 

maww»aw» 



4.1 Foz isu3^t:ioa of £h© Equations of Hatton 

In en inerei&l coordinate oystesa, the ©ofcton of tbs 

vehicle am bo described by equating its acceleration to 

its grevitsfeicnal specific force. Since it is decided to 

•» 

describe the motion in the coordinate set v-mish rotates x^ieh 
tha maon, the additional apparent accelerations which appear 
in tbs coordinate traaaforsaafcion ausfc be included. The total 
acceleration consists of tfea saas of tha following terns: 

1* The acceleration of the coordinate system in 
inertial space, This is considered aero in. the restricted 
tforee-'body probless* In the physical problem in cislmar space, 
it corresponds to the accelssrat&aa of the aertfe-sscch oy&te® in 
its orbit iisotsid the mm, including the perturbation effects 
of other planets. Its only significant component is ths cen- 
tripetal acceleration of the narth which sees© to fee the ssost 
appreciable effect neglected in confining tha investigation to 
eh© restricted threa-bedy problem. If later it appears desirable 
to consider this ter®, ie can be included as a periodic forcing 
function in ths equations ©£ motion, representing the differer.cQ 
between the specific gravitational fort© of the sun ©a the . 

Uk 



vehicle &nd that force on the barycesatsr of the e&rth-saeon 

system. 

2. The tangential acceleration. This is also 
considered ho fee sero since the moon 1© assumed to be in a 
uniform circular orbit about the earth. If the effect of 
the eccentricity of the t&oon*8 orbit were to bo considered;, 
this ter© would have to be included. 

3. The centripetal acceleration. This la treated 
as a centrifugal specific force and grouped with the gravita- 
tional forces. 

4. The relative' acceleration. In the (js, y) rota- 
ting coordinate system* its components are k and y* which are 
included in the equations of motion. 

5. The Coriolis aseeleratioa. Xts magnitude Is 
vjTTs ~2 

Z&vjfc «J* y , and its direction is defined by the vector 
cross product of the system angular velocity and the relative 
velocity of the vehicle. Sts 2 $; and y coaspca snfcs are -2Jty and 
4’2&fe respectively. 
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Fig. 4.1. Representation of the Rotating Coordinate Set. 



As sfesswa in Fig. 4.1, the origin of the coordinate 
aystesi is chosen at an equilateral trieagle point, as 
Trojan point, since t?e are primarily concerned with the 
dynamics c£ a vehicle in the vicinity of this point. 



The ss: and y exponents of the vector equation of 
motion are then 

* - 2n * * £ lx + £ 2 k + f cx 



y + 2fl* - f ly + f 2y + £ ey 
A 

Substituting gross equations 3.1, 3.2, and 3.3 and resolving 

■» 

x ~ 2Clp to cob y 4» cos £ 4 s 6{®j4®^) r cos 

T* ST 

r l 2 

y + 2ffi „ ('IGCaj^eSg) ain ^4. ^CfSj+Og) bJj* (5 + GfeMa.) r ni» 
r ._ r 



1 



2 



To further simplify Che equation®, time will be normalised 
with respect to the period of rotation of the earfeti-sason system 
so that 

j.2 a G(ai|j4Qg) j, «2 

’3) 



where t ’ i3 now real ci«no 



with fcho dot denoting differentiation with respect to 

dimensionless time. 



In dimensionless time, 
becoEct© 



x * 2 y 



n 



the equations of motion 



cos /A* r cos ^ {4.1} 



1 



y *4* 2 x ® - — 4*~» sin >j' » sin/ 5 *!- r sin ^ 

r l s 2 <^2) 

Referring to Fig. 4.1, the quantities x^ 9 r.,, and the 
trigonometric functions of ^ ^ , and y sa*^' be expressed 
in terms of k and y as follows: 



r? ® Cy 4* V “) 4* (s *s* -4) « y* *s*^3y f / 4* k + 1 



ilfl 



/ o 2 



2 .C 



2 

A 2 



- 2 



2 



r 2 «s (y -f -•) * ( -- - s) « -f Ay + x - x + 1 
2 - 2 - 2 



o f« 



r cos T B3 r: ■$* /) - - ■ K i -4 

2 - 2 



t 



y * 



A 



COS y 



sos 



Sin fj' 



n + -i 



£3 



r 



1 



(s 4* %) (y 2 »’« Ay + xA x -I- i) 



P a J^JLi 









2? 



A 

V J.W 
> * 



2g) (y 2 + Ay 4* ?£^*» x + 1) " 



as 



i w (y 4> V <|} (y 4, 4* Ay 4’ X 2 + X + 1) * 



in ^ *-• 



y 



7 

* ^ 



' 2 



r _ 3 

(y (y 2 + V 3y 4- :; 2 - x + 1) 



Substituting in equations 4.1 and 4.2 the equations 
of motion become 

j*t-2y *= x-v'' 2 -G “ (ji-f'i) (y 2 + 3y+3j 2 -f:;+l)/2 (y 2 4-Jy4*x 2 ~x+l) ^ 



<*.3} 

y*+ 2x « y4- ~ - (y+ ~) (y 2 *i* oy+-x 2 +:rf*l) /^-' : (y+ -) (y 2 4$Ty+sc 2 -»HL) 

2 2 2 

(4.4) 

4.2 Integration of the Equations of Ilotion 

The only known integral of these equations was discovered 
by Jacobi. It has since been proved that, in this coordinate 
system, no other algebraic integrals exist. 

Let u «= “Cx+'a(y 2 +Jy+-x 2 X'41)4* 4 (y 2 '^y fc ^* 2 +x^l) 

_ J. 

*K(y 2 *^y+x 2 -»KL) " 

(4.5) 

Then ~ «= (y 2 4%3y-f*x 2 4*:^l) 

dx -3 

- e (x~ h ) (y 2 + ^y-i-x 2 -:ri*l ) 2 

(4.6) 

“ *= v 4.1 *- 71 (yf — ) (y 2 4"'Cy<-x 2 H-:c4-l) Z^ 2 

ay 2 2 3 , 2 /- 2 : 3 / 2 

~ G (yb — ) (y v Jy4*x -:c4*l) ^ ^ 

2 



and the equations of motion are 



2 y 



du 

dx 



(4.7) 



(4.3) 
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<4.9) 



y+ 2 sk « — 



By 

Multiplying eqaatioas 4*8 *ftd 4*9 ky 2as and 2y respectively 



<** * 

2 i s » 2 



• « 

22 S3 2 % 



M 

Bs 



2 yy * 2yK 




Adding equations 4*10 end 4*11 



*• * *>* 9 • 

Ssrjs ■;* 2 y | » 2x 




+ 2y 




(4.10) 
(4* IX) 



«* d& . **« dy _ « Su ds , « Bu dy 
d£ - , dfc dfi - dy dt 



Xntegr&ting 



j * "v- ^ ® 2 diE *{- 



Su 



ds: 



• By 



)dy) » 2u v 3 



e •> -2o + i 2 +y 2 



The Jacobi integral is 

Q ra gejx «■- ) °2 Cy^3^a^arf»X ) ^ 

“2® (3^4*^3yi«^«*SB>l p* *•!* «5» y^ » ecsaofcsnt 

(4*12) 

4*3 Vectorial Fepranca&fcafcioa of the fiauatlaag of Mat ion 

— n *m> 9 » »Bp*«we irr ntr~ ir~i~ mi — 1 ~- r-r - ~ -"'ir— riTr~~~Tr — r r r n — r~~ 'i r tt— -; i — t 1 t* — r ttt rrir^ri 'T' r~ n n nn mnnTTinflrnri i -Tir rrnr nr rrri mfrrurrr 



It akculd be noted h/sre that t ho sufesequesst 
oqmttffiw apply to osaly one of the two Trojan points* For 
the other Trojan pointy & syima&vy about the as-axle sppl iee s 
but the oqpM'l&m and direct ions rosy not be applied literally 
without eoasidoffdfcicsx of this ©yasa-atry 

IS the asotion oi the vehicle Is represented in the 
eofsplox plane where 
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s ® s 4- iy 



equations 4*3 cad 4*4 sass bo ■sjrittsn as 



2S 



2 y » f _ 



(4,13) 



iy 4* 2 is ® if, 



y 



Adding equations 4*13 «ssd 4,14 



(4*14) 



a 4* 2 is » f R <** i fy 



(4,15) 



At £isefc thought# it would sees that* for a stable 

Trojan point orbit, fifes resultant of £„ and f y should be 

directed fas&ard the origin* TI&b is not tbs o©se# Is will 

be seen later filiate in fifeo vicinity or tlie trojan point# ^ 

sad £ have fife© eaos sign ae x and y respectively* aatfeer, 

y 

it Is tbs Coriolis acceleration tthlch provides firs average 

relative acceleration, s’ toward the origin if the orbit is 

to be stable? The vectorial representation of Eq* ^*13 is 

aSscssn to Fig* 4*2* This indicates that tbs orbit, sutt have 

a gross stotlon elocttwise about the origin for stability* 
in- h 

8 




Vi? 






.*»»»■«**$£> ^ 



Fig* 4*2. Vecfiorial Representation o£ the Bqu&ttas of Motion* 
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4*4 Linearisation of the Equations of Motion 



As a first order approximation of the behavior of a 
vehicle in the vicinity of a libration point, the equations 
of motion are linearised and fcha solution of the linearised 
equations is examined. The procedure used by Moulton (3) 
is essentially to linearise the right hand sides of equations 
4.1 and 4.2 by Taylor series and solve the resulting linearised 
equations,. Whittaker (4) uses a different approach, working 
with the equations of taction in the fora of four first order 
differential equations, with the sssze results. 



The first approach will be extended here to obtain a 
more complete linear model of the system, including an integral 
of the equations equivalent to the Jacobian. This will later 
be used for a more complete comparison of the actual and 
linearised systems. 



Linearising the differential equations is equivalent to 
retaining second order terns in the Taylor series repre- 
sentation of the Jacobian function* Eq. 4.5. Letting primed 
symbols represent the quantities of the linearised system 
Eq. 4.S can be written as 



u’ <s,y) 2 u ( 0 , 0 ) -s-^a (0,0)x-5-|S (0,0) yH'--§ (0,0)s S + 



(4.16) 



d 2 u 



^2 
O U 



(0,0) (0,0)3/ 

cSKaj $ 3 / 



Substituting x « j » 0 in equations 4.5, 4.6 and 4.7 
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V 




U { 0 , 0 ) \ S «s 



(C,0)» -s + ij - % e « 0 

-2a (o.o) - 3 -^5 - « r l =. o 

ay 2 *2 2 

Differentiating equations 4.6 sad 4.7 

i+3^(^) 2 <y 2 ^yfx 2 +»»*l) ^ -^(y 2 ->€y^K 2 -5*5£4*l) % 

8x 2 

+3e <x-%) 2 <y 2 *^y*x 2 -x* 1) -e (y 2 + v 3y*x 2 -xl-l) ' 2 



<4,17} 



.2 



« 1*3?} (y*5) * (y 2 * /Syfx^f xfrl K -? (y 2 * vSytaA^l ) ? ' 2 
dy^ 2 

«^(^^) 2 (yVW -*uf?2« (y 2 */3>K 2 ^l)“ 3 /4 

<4,18) 



® 3^(:e*?2) <y*^|) (y 2 *{^^K 2 4«*l) 3 6 <x~Jj) <y/l) <y 2 *Syl*^-r*X)'2 

Sxdy - 2 

(4.19) 

Substituting x » y » 0 in equations 4, 17, 4.18 end 4,19 
2 

(0,0) ® 1 + ~ ^ ^ e-€ B 2 

dx 2 4 4 4 

2 

<0, 0 ) *» X 4* ^ *{* ca -2 

3y2 4 4 4 



>2 
o u 



3<i 



dx8y 



4 



'£l <L m Q - 6) • 2.5 <l-2«) 

4 4 4 



Substituting these quantities in equation 4. 16, the Jacobian 
function for the linearised model Is 
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V 















8 3 3 2 , 

® ^ ~ s *1* 

d a 



(X * 2®) xy + ^ y^ 
4 8 



(4,20) 



Di££e2 , ©Qti£i;iE2gj > the linearised equations of mot lea or© 
obtained 



s - 2 y 



cha* 3 * 3\/S 

5s 4-4 



(X - 2«) y 



(4.21) 



y -r 2 X a — * jil (1-26) & 4* ~ y 
Sy 4 a 



(4.22) 



cad fc&a Integral of the aquations* equivalent to the Jecobiea* 
is 



e* 



* 2 
» 22 



*?» 



* 2 

y 



3 





<X«2$) xy « 



(4.23) 
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Gwmm s 



to Bt&mtm m to mu&iim gj msic® 



5 . 1 T he Stabili ty &r&£&gio n 

wBW»WWMfr»iniw. * f» i"4 ,i ii>ki. w7i_uiiiM>L > <*^*ffia»>-a<Kaccwr>ciwi>iii>i irwiatwi m*?’ 

In cpas-aSes- £osta» ©qusticBS 4*21 

and 4*22 ess. bs written as 

(p 2 - — -) s « j ap *5* ~ (1-2«)| y <* o 
4 , U. 4 ■"' 

|2p ~ {2-2 e )j ss *4* {p 2 *» ^) y «■ 0 

4 _J 4 



The ehaK&etcsflBtlc equation o£ this eytten Is 



<P‘ 



2 



*) (p 2 - ~) 4- Tap 4 — <1-2 e}] [sp - — (l-2G)j« 0 

y 4 *— 4 - J ^ 4 



oaf 



P 4 4 p 2 4 — s (l«e) « 0 
& 



The toote of £h;s shetaoterloslc equation a^a 

era c * £j 4 'j $$ *» ««£ § 

- 4 

S£ the second tom io seal* ’She cfcaffec&et&stic equation 
yields t®» paiirs of pure imaginary toots* and the solution 
consists of pute stecsoido* Scow/tssr^ if the second tern is 
tsMSlnasy* the toofts hoRw pos&t&re s?eal p&srtSa and the cm~ 
plitude of the solution iasresses eRponeotielly with fcise* 
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Ifoa stability criterion is s$vm 



OK 



% ~ ~£ e (l~z)>0 

«* 




GU 

In tha ©srth~®ooa eystssa* -« is 

ffl X 

Che stability criterion is esL 



1 

**»**« n 

81,3 



anti 



$«•£ th e Sgggig!|S«Liafi Libratioa Point© 



Stoultoe m®Q «s similar dsv&tasasafe to shew that sh© 
straight- line libratsioia points qsq unstable, Ano&h <ar presen- 
tation of this approach Is ssade by Buefca&a (3). 



If tbs equations of motion ere linearised about a straight- 
lisa© liferatlem point 9 tbs ©hssractioristiei equation of the result- 
lag ©ystesa yields two toeginary roots and two real ro ot8 a 
one ©f which is always positive* A part of the solution 
must than ineireas© e^enencially with time* The aystea 
ie tbue coasld©r®d was tab lei a assail disturbance in dispIessssEat 
frcsa the libretion pains will increase with tiss. 



3 >3 fffca ggols a F&te sS Solution 

*n rm urmn w rai qe*5iW^ , 7xcr3^g»g^re mm reKUMiu^stu* 



2 



ta_ 

is snail eacr?ji?^s., as it is la cfee 
l 



IS tfhst m@® x&k'lo 
©&rfch**o«& ©y etcmj felt® Xtos&jrised eqiu&lo&s of s»fci©& yiald 
a solution of tks £&m 



s « sto p^fe *5* go© Pjfe 4 * G- sin p^fe ❖ cos p^fe 

(5.1) 

y * Au) sis pgfi + Bg oos pgfc + Cg 3te pg£ + D« coa p«6 

(5.2) 

■» 

ft n 3?>r 'eras 1 m **•? * 0 T 

P$ ® % 4- % VI - 27« (!-«) 

U> 

pf * % - %\ r rr~2?€ (3^«) 

end fete? of fcfee eight ecssfeanfes of Issfeegsafeioa ase independent* 
ffoo r «X*fc ioneMpe moag the constants csss be dafesfemtoed by sub« 
&fei£u£tog the solution b&sk into tba Xiaeasr diffaveafiial ©qua*' 
felons* Solving ffosr she eons tense of Eq* 5*2 to tee of tbs 
constants of Eq« 3*1 




- 2d. 
-All. 




(5.3) 
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is A) 



' s 1 ’ I fd l -?A) 






2 



^ 9 

r*A 4* 

P X 4 4 



C rt w 



4^3 (l- 2 e) €3 •• 2 p g B 1 

• XM- **- XWO'U^J*. 



+ 9 



4 




C3* 




r'^'r '<-' m — — *"*csrs*«.'V 




( 5 , 3 ) 



(5,6) 



The coasts- M:* A- , B-*, 0 , » Is avre than datetnlnad by 

, 

the initial ess^litieno of displacer jesac *ad velocity, ss (0) 
3ttd >' C®l * 



* p Ar ifoy n 

y ^ 



5,4 BliscinetSc?* of. €te» Mi'xilk o£ O&gj L i'Ae t ioa 

Ji <M# WWIIWl »MB ~ **■,■ **\ - g - ,.>ujirw - iwlWvi- - - •-'»' .1 ■ ~M~. . 'UC-...-HO. ***,' »« _ 

to better ujiiesgtead f.i*» •?.•. irure -of tfe© solution that 
has bce-.» obtain -? fro* ths Xfnaar&Q&d equations, it would 
be deeirable te dstesisdbae hors fcb© nolotion 1* divided bow,, 
the ®Kr itodae of vibration r : a ftufrfcioa of initial eo»d£t'.ic->o . 
Tfeo ccKStcuot:© ef &%uet£cna 5.1 *?.d S.2 can be detemia^sd *? 
a fuaetier of initial csmdit&oao, :ia feera* of tin natural 
£xeq&m>.o&t.9 p-i and p«. This to d«j« la /upp«fedis A using 
snaMsiesl ^oodCrnta i:er the cuvtti-rxo.i oysfcesi* but tfca 
physical prcc^cr by which tha sysfc.jy. is excited to differ* ni 
proporh&asc of the two «od®o of vibration ic not readily 
app*rc3\fc. 



To obtain further insight along this Xisi<%tfe determine 
what specification of initial conditions will excite only 
erne Esods of vibration, say the p- node. For simplicity wa 
arbitrarily specify the initial displacement as 

s ( 0 ) w 0 ( 5 . 7 ) 

y (0) » 0 ( 5 . 8 ) 



We are, of course, concerned with crash smaller displacements 
but the results can be scaled down. If ws specify that the 
initial velocity is chosen to eliminate the Pg sede, then 
equations 5.1 and 5.2 become 

x *» A, sin p^t + B 1 coo p^t 

4- Bg cos Pjt 



y a A 2 82Ln P^t 



Substituting initial conditions 

X (0) a, B 1 ® 0 
y (0) « Bg 03 1 

x (0) <=» p^A-j 

y (o) ® p x a 2 



Also from Eq. 5.3 




•» 



3 

4 ’I 



3 (1-2S) A x 





4 
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b, - 

2 ju | 

Pi • z 



Since Bg 1 



chan A l - JL (p 2 + 9 , 
2P l 1 ^ 



and 



* 2 “ * ip | 0 - 2 ® ' 



Thus the initial velocities mist be specified as 



A (0) 


® h CP^ *&* 


— ) 


(5.9) 


1 


a 


y (o) 


-- 


(1»2C) 


(5.1./) 



And the ouboequsnt ssotic® of the particle is jpesified a© 

(5.11) 



n 2 ,9 

» P 1 4 ©in p^t 



2Pj 



CS 



y 63 



I_1U1=2S1 ain 

2p x 



4* cos p-^fc (5.12) 



Solving for sin P|£ and cos p^fe in t©ms of js and y s 
can §h*n eliminate tics* by eh© equation 

j[s£n 2 Pl t] (x,y) *J* cos^ p«t (x,y) <■ 



After simplification., the equation of the "orbit” about 
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* 



the Trojan point is found to be 



-,2 lllxi-lfi 

<1 J. fl \ 




2 



2<p? 4 9) 
L 4 



(5.13) 



This is the equaticn of an ellipse with its center at the 
origin. It shoulc be noted that this is not a Kepler ian 
elliptical orbit. Angular momentum is not conserved in 
this coordinate system, and there is no central force 
field analogous f.o that of a Kepler ian orbit in inertial 
space. 

Although th? motion of the vehicle about the Trojan 
point described to the relative coordinate system may not 
be correctly ter ted an orbit, the term will be retained 
for convent ©ice. 

If wo let v be the magnitude of the velocity of the 
vehicle in the rulative coordinate system, i.e. 



Then this velocity can be found at any point on the orbit 
from Eq. 4.19 
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C Is a constant which can be calculated from initial con- 
ditions 

G' «« v 2 (o) - 3 - ~ 

b. 

Fro© equations 3.9 end 5d0 



v 2 (O) » \ <p 2 + + 27 <l-2«) 2 

64 



«nr* at any point on the orbit 



V 2 ' * % (t>? * 9) 2 + — <i-2€) 
1 64 ’ 



4 



~ 4* - 3S 2 4* -/g 

4 4 2 v 



(i-26) xy 






(5.14) 



Considering an arbitrary point on the orbit to be the 
initial point., the initial velocity necessary to excite only 
the p.| mode is determined. Its mags aitude is obtained fross 
Eq. 5.14* and its direction is opacified by the tangent to 
the orbit computed fron Eq» 5.13. 

If we nmi divide all displacement and velocity varia- 
bles la equations 5.7 through 5.14 by an arbitrary sealing 
factor* we obtain a new orbit with the sase center, shape 
and orientation., which is also a solution of the linearised 
equations with initial conditions specified to eliminate 
the pg mode. For a continuous variation of scaling factor. 
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wo tm infinite ismiiy ©£ ellipses eovesr log the »~y pte, 
25(0) «* ) 0 (3, IS) 



^(0) es k 


<0 , jLfe ) 


K(0) » ■& (p‘| 4* ~) 

2 A 4 


0.17) 


y<0) e» ^S\«r (1»2«) 
8 


<5. IS) 


52 cb JL. (p? 4* ~) Sir* p^t 
2p^ - 4 L 


<5.191 



y s» •£& ^3 (1-2*2) sin p-5 4 * k so 3 
8^1 

P -n r 2 

2, I 3(3 q-2c)~ 



y 4 



2<P? 4- 1) j 

~ „ 4 



*S* 



4p 2 ^ 



16 (pf 4* ™) 

L. A . 4 



9 

\ m , 

h *» k 



<5,20) 
2 

(S,21) 



2 



2 * *■"* 2 9k 2 , 3 2^2 + 



V” tsv ^~>(p? *|» S) •§» (l*»2 e ) 

4 1 .. 4 . 64 



4 4 



~> f3<l-2«)ss+ - 
2 4 




(5,22) 



Given «&y point in the x, y pi gas 30 fcha initial point* 



33 







% 





















we c @r calculate k by substitution in Eq. 5.21. Tve can 
then determine the required initial velocity vector to 
excite the mode only from equations 5.17 and 5.18 and 
the resulting orbit is specified by equations 5.19 through 
5.22. 

If we wish to eliminate the p-^ mode instead of the 
Pg mode* the procedure is the same* and the resulting 
equations are identical with equations 5.7 through 5.23 



We now hav© the mathematical equipment for answering 
the question, "Given an arbitrary initial position, what 
vector initial velocity is required to excite only sod© 



sent at ion for answering the question "Given an arbitrary 
initial position and an arbitrary initial velocity, what 
does the orbit look like?" The temptation is to eubsti« 
tute mashers into the equations of Appendix A. However, 
we can do it sssdL-grsphieslly by a two stage trial and error 
method, and perhaps gain & better physical incigfet of the 
problem. 

X*q 6 'tlsa equations of the orbit be written as 




5.5 Determination of the Orbit from Specified Initial 
Conditions 
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X ■ + Xj 



< 5 . 23 ) 



* m (5.24) 

whore ■ A^ sin pjt + cos p, t 

y-^ » Aj sin p^t 4* B 2 cos p 2 t 

k 2 ** ®i 8in P2 C * D 1 008 p 2 fc 

y £ " C 2 sin p 2 t + Dg cos p 2 t 

To determine the orbit resulting from specified initial 
conditions, pick a point £acj(0), y^CO^TJ as a trial. This 
determines an ellipse of the family of Eq. 5.21 and specifies 
the required x^(0) and ^(0). It also specifies the point 
£^ 2 (0 ) , y 2 (0)] by equations 5.23 and 5.24. In turn, this 
point also specifies an ellipse of the family of Eq. 5.21 
with p 2 substituted for pp and gives the required ^(O) and 
^2<0). If the trial is correct, the vector addition of v^(0) 
and v 2 (0) gives the specified initial condition v(0), and 
the orbit can be constructed by the addition of a vector 
generating the p 2 ellipse at a relative angular velocity of 

& 

Pi * 

This method is Illustrated in Fig. 5-1 which is not 
drawn to scale. It io included only for possible clarifica- 
tion of the above explanation. 
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Desirable Ranges of __Xsa.it ial Conditions 



From t^e preceding analysis, it; appears fchet, given sin 
unavoidable displacessent error, we will obtain the ''smoothest 1 ’ 
orbit if the initial velocity has a magnitude in a range such 
that the Jaeobiax^conatant of the vehicle is approximately the 

' "" * V . 

aaaaeTa* that Of a body at rest at the Trojan point- It can 
be seen by substitution in Eq. 4.19 that a body vibrating in 
the t 20 ^ only has a relative energy higher than that of a 
body at rest at the Trojan point, while a body vibrating in 
the p£ mode onl?’ has a relative energy lower than that of a 
body at rest at the Trojan pop£ht.‘. 

More izsportanf, we wish the initial velocity to be in a 
direction compatible with a clockwise orbit about ths Trojan 
point (using the notation of Fig. 4.1). 

The method o£ the preceding section shows that, for ini- 
tial velocities in the "correct 9 ' direction, we tend to get 
contributions from the two modes which are additive, and the 
outer bound of the resulting orbital motion is not much different 
from the initial displacement . On the other hand, if the ini- 
tial velocity is in the opposite direction, we tend to get 
large contributions of the two modes which subtract at the ini- 
tial point, but add at the outer bound of the orbit. 

This verifies the impression we get from section 4.3. 

Since the gravitational and centrifugal force field is, in general, 
directed away from the Trojan point, the Coriolis acceleration 



1 
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cs«st provide the net acceleration toward the Trojan point 
for a stable orbit. For this to occur, the gross motion of 
the vehicle must be clockwise as can be seen in Fig. 4-2. 

Fig. 5~2, whicl is also not drawn to scale, is a sketch 
of two attempted in, actions of o vehicle into a Trojan point 
orbit. In each case the approach is from the negative x 
direction, and the vehicle has been provided with approximately 
the correct relative energy. The closest point of approach 
to the Trojan point is considered the initial point for 
generating tbu orbit. 





Tjojan Point Orbit 
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Under these conditions* a miss to the left of the 
Trojan point results in an orbit whose maximum displacement 
from the Trojan point is much smaller than that of the orbit 
resulting from a suss to the right of the Trojan point... In 
an actual attempted injection of a vehicle into this type of 
orbit it appears we should evaluate the expected angular 
error in guidance ^ and then aim to the side which will insure 
a resulting orbit of the first type* 
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CSfiFSWL 6 

«<|WI 1 . H T01W..I I n 



GC^IPARXSOH OF THE AC 7 CTIAL BQUAfflO^S OF MOTION 

mm THE LIK&RI2ED EQUATIC&S 
* 

W© will define, T, /.he residual force, so be the 
vector difference between the actual fesee exerted on a 
vehicle by the cabined grs vitafcional and centrifuga? force, 
and the fictitious force w.Vidh would be exerted on the 
vehicle by the linear iced node I of the system. 

l w a 1 ~ T «» gjad « - grad u 8 » grad u’^ 



where 



* 

u 



S3 U~u' 



and the quantities u and u s are defined by equations 4.5 

end 4.6. 




> 1 



Figure 6.1. Effect of the Residual Force on the Acceleration 
of tlv ; Vehicle 
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The sre/iidusl force directly represents the effect of 
all higher order terms of the equations of motion which 
were aegXscted iu ehair linearisation. As suown in Fig. 6* 1, 
if the vehicle Ip assumed fco follow the solution of the 
linearised equations a fifes approximation, then the reel- 
dual force vector shows she tendency of the actual system 
to cause the vehicle to deviate from this path. 

It would be convenient if the residual force field 

were found to have sees s significant pattern. For exsgspXe, 

if a solution of the linearised equation© was formed where 

the residual force was always inward cow&raths Trojan point 

•* 

this would be an indication of the existence of a stable 
orbit somsf where within this region* 

Several points of the residual Jacobian function, u t 
were calculated by the sasehodo shcgm its Appendix B» The 
results oi these calculations are eliowa in Fig. S.l* The 
residual force is the gradient of the function u‘ t i.e., 
it sets along the direction of greatest rate of increase c£ 
n with distance, and its mgnifcvde is equal to that \rat^ of 
increase. The plot of Fig. Q«X is neither sufficient iy 
complete nor accurate Plough for a numerical determination 
of the residual force. Xn any case., it is more easily ealeu 
lated directly by components with £he methods of'Appaadix B. 
Instead 'die plot ms meant to give a rough indication of the 
direction and magnitude of the residual force as a function 
©£ displacement from eh© Trojan point. 
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tio significant pattern of residual force distribution 
was found. If we trace out an orbit about the Trojan point, 
we find that the residual force flunctuates in magnitude 
and direction with no readily observable net result. We 
also see that the closer we are to the Trojan point, the 
smaller is tbe residual force. 

So far, we have found nothing to indicate that our 
linearised model is not a good approximat ion of the actual 
sy 8 tea. The deviations show no significant patterns, and 
they seam to decrease regularly as we get closer to the 
point about which the equations were linearized. However, 
we are still not guaranteed that an orbit about the Trojan 
point is absolutely stable for "small disturbances", nor do 
we hove any information about how small the disturbances 
must be to obtain "partial stability," i.e., keeping a 
vehicle within specified limits of displacement from the 
Trojan point for a specified period of time. 
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CHAFEEa 1 




GmmmTim of m cam 



Usually when we seek m analytical solution of a systea* 
of soa~ linger differential equations, we fed that wa esn 
always fall bae-k ©a numerical method, in nil else fails* ¥e 
would like to have a solution in closed f« or, in this esse,. 



ea analytical upper bound ©a the orbit, but, ia general, this 
is © convenience which saves us the tSsa© and expanse of resorting 



to a digital eampu&er* 



Here, there la setae. doubt about Aether 



we can obtain a satisfactory answer to the stability question 
by numerical, methods. 



The obvious method is saerdy to Integrate the aquatics?© 

©f motion numerically over ae long a fcisa interval as is required,, 
and thus obtain a numerical representation of the position of 
the vehicle throughout the ©peed fled time interval* The limita- 
tions on this method fall into two desses® 



the first limitation 
physical constants of the 



is our inaccurate knowledge of the 
syst«2sa» Although it Is believed that 



our estlmtes of the ©artfe-mooa mass ratio .way fee in error by 
a® saieh as 0*3 per cent* ©us present Xiaor rocket ©uporiamfe® 
should significantly advance our knowledge in this area* In 
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asy cose, this i* »©fc a basic lissitafsioa of the isxves&igafcio a 
of stability* If we can obtain she answers to th& stability 

SR* 

questions based on a model of the £ssrth«ax}®a systesa where 

is sssactly 0*01213, we will hsore taken & long step toward 
obtaining answers for the actual systees* If later, with more 
precise knowledge of she mass ratio, we feel that She results 
do not apply to the actual system, the procedures for compu- 
tation will still apply. 

The major limitation on determining a precision orbit is 
the errors inherent in the method of computation. The affect 
of the choice of a particular method of computation is dis- 
cussed by Baker and Others (6) . Basically, these error© asey 
be classified a© truncation error and roimdiug error, Trunca- 
tion error results from the use of a finite interval over which 
the value of ths integrand is apprejsfcgaeed* It can be reduced 
by decreasing the else of the interval, end thus more closely 
appreciating the value of the integrand over the range of 
integration* Rcimding error is due to the limitation of compu- 
tation to a finite maabsr of digits, sad She error in each step 
of She ecsBpuS#l*.ic*a la propagated throughout the satire problem, 
Ths resulting rounding error then increases with the number of 
steps » 

These errors cm bo reduced by intelligent selection of 
methods of eoEspu&atioa* For example, for an orbit of this type, 
it is felt that computation of perturbations from a referee© 
orbit would be more efficient than, direst integration of the 
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equation*.. H ©waver* the cc^ijuttion of truncation «m« aad 
rounding error still can not be eliminated, For any precision 
orbit eoaputatioiij there is a practical limit So she deration 
of the orbit which we eaa compute without a loss of validity 
from the cumulative error, Thus* for a general luaar probe., 
regardless of refinements the method of ce&rantetian, ws 
csa specify where it will be on its first pass near the moon* 
we can estimate where it will be a few smtfcs later* but we 
have no idea where it will be several years later. 

Intuitively* it is felt that we emi do better tfes» this 
with &. periodic orbit such m ours. It seems that* If the 
orbit is periodic* we can confute a precision orbit to comp loti-', 
one revolution about the Trojem point* sad the orbit S.e stable 
if: it passes closer to the Trojen point the second time around , 
But wo feevo neglected the motion of the vehicle in two ©odes 
with inoossasnsurebla periods* A hypothetical vehicle moving 
in the linearised system will ov will not meet the ©bov© cri- 
terion according Zq its initial condition®* hut W3 knm? that 
the orbit of this vehicle is stable. It raetaias within upper 
and lower bouada determined by tlm initial, conditions* and will 
pass arbitrarily close to say point in this region if ws wait 
long enough, 

the ue%£ thought is that we should then compare tfe® computed 
precision orbit with the solution of the linearised equations « 

But again this comparison gives no conclusive evidence of sta- 
bility* At any comparison point* fc&o motion of the vehicle in 
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the actual syslar* will hm/a altered site proportions of rels- 
tive energy between the two taodes. Wo my find that the ac- 
tual orbit is closer to the Trojan point, fcsufc if we ccaapute a 
new linaarised orbit using an initial conditions the actual 
displacement and velocity at the comparison point, we may find 
that the predicted upper bound of subsequent displacement fraa 
the Trojan point is greater than that predicted by the first 
linearised orbit. 



We o an, of course calculate the predicted upper bound and 
continue, watching this upper bound for an increase or decrease ^ 
But now we are back to the problem of: not being able to watch 
it long enough, because of the cumulative error of the numerical 
calculations* It is this cumulative error that nay ©ask the 
stability affect that we wish to determine, especially If the 
instability of the orbit is ssa&XX. Me ®m working in a regie®*, 
of ciG lunar opaca roughly bounded by two ellipses defined by 
fch® solution, of the linearised equations. Sine© we see no sig- 
nificant pattern in the distribution of the residual forces (or 
non- linear effects) in this region, we should cover the entire, 
region rather uniformly in ©ur calculations • Otherwise, ws 
cannot foe sure that the navt £cs» revolutions of tbs vehicle 
about the Trejan point will not pass through points where the 
non* linear effects ere entirely different* 

She length of time necessary for m orbit to pass within 
a specified distance, <f , from every point, in its permitted region 
varies with tbs distribution of the motion between fcfee two nodes 
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of periodic motion, If Shs motion of She vehicle is primarily 
confined to on® sods wish the of;h«r ©£ amplitude no greater 
tbso <f , then the necessary fcima is %x o longer than the period, 
o£ the pr£sm*y s&ode* If the amplitude of the secondary mode 
is large 3 the tlffia may bo prohibit lira due to emulative error « 



This leads os to consider invest ig&fcic® of the special 
cases where* based on the linearised mode 1 of the ayefcsa* the 
initial conditions are such as to eacits only one sods of 
vibration- The trouble with this is that we should expect 
Che non-linearities of the actual system to cause transfer 
of energy Co She other roods * with a cumulative effect until we 
are back with die general ease again* 



Xfc seesaa possible that wa may obtain useful infomatien 
by constraining the orbit So the simple ellipse,, mid inte- 
grating the effect of the residual force over ones resolution 
about the Trojas point. Intuitively., it is fait that a net 
inward fores would be an indication of stability while a net 
outward force would indicate instability* X£ this ware done 
for orbits of various sisss, we would then have an indication 



of stability versus average distance fro© the Trojan point * 

This investigation has disclosed no mathematical juatificatlaa 
for asserting the iseaniagfulness of this "stability indication* 
However it is felt that the discovery of & valid stability 
criterion will require calculations of this type- 
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CHAPTER. 8 



LINEARIZATION OF THE EQUATION OF MOTION 
AT POINTS DISPLACED FROM TOT, TROJAN POINT 



8.1 Linearization o£ the Equations of Motion at the Position 
of the Vehicle 

In a digital computation of the orbit of a vehicle about 
a Trojan point. It seems most: efficient to use a reference 
orbit and compute deviations from this orbit. For a reference 
orbit, the first choice that comes to mind is the solution of 
the linearized equations as developed in Chapter 5. However, 
for a small segment of the orbit, we should get a better fit- 
ting reference orbit by linearizing the equations of motion 
at a point on the orbit, rather than at the Trojan point. 

Linearizing equations 4.3 and 4.4 at the point (a,b) 
rather than at the point (0,0), the Taylor expansions are 



5u 



2 9 

— (a.b) (a,b) (x-a) + (a,b)(y-b) 

dx ax” dxdy 

( 8 . 1 ) 



2 2 

cs (s,b) 4 * — H (a,b) (x-g) ~~M(a,b) (y-b) 

by by bxby oy^' 



( 8 . 2 ) 



For convenience in writing the equations, let Rj and Rg 
be the distances from the point (a,b) to the earth and moon 
respectively- Then frcm Figure 4.X 

R? » (b *S”^) 4* (a 4* %) 2 =* b 2 4*'f3b 4" a 2 4* a 4* 1 (8.3) 

1 2 

jt 2 2 

R 2 » (b 4* -) 4* (a - h) » b 2 +«Sb 4* a 2 - a 4* 1 (8.4) 

2 2 

Equations 8.1 and 3.2 can be written as 



du 

dx 



° k^ 4* kg (x~a) 4* kg (y-b) 



du 

dx 



k^ 4- k 3 (x-a) + k 5 (y-b) 



(8.5) 

( 8 . 6 ) 



where, from equations 4.17 through 4.19, with substitution 
of equations 8.3 and 8.4 



k, « — (a,b) = -A (aHs) &T 3 -€(a-*i) Rp 3 

dx 



(8.7) 



k 



. a 2 u 



.3 -3 ^ . 2 -5 2-5 



2 ~^2 ^ b ) ~ + 3 ft (art*?*) R 1 + 3e(a-^) Rg 



( 8 . 8 ) 



kq = — (a,b) = 
dxdy 



r *-5 jt ~ 5 

3A(aH)(bK“) R-, 4* 36(a-2s)(W* f3 ) R 9 
2 *■ 2 ^ 



(8.9) 
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Sm <fh /«a “3 7 o *»3 

88 «5 4* ^ (b 4* «~) R-« (b 4 --) Rn CS«IO) 

Sy . 2 2 L . 2 * 



k e 



i2 "3 "3 V"V 1 »s /o 

* «-§ <s,b) *» 1-^Rj, - *SU ♦ 3*(f> * ~) R x 3 r3«(fc+ f) 






Ur 



’5 



\W» 



£a litsesrised £ 02 ©* equations 4*3 and 4.4 cm. bo ^srifctee 

an 



»* • 2 y » (&£ * akg - fek^) 4 k^s 4» k^y (3*12) 
y 4 2 k » (k^ ~ ekg - bkg) 4 - k^s 4 - k^y ( 8 * 13 ) 

Tkesa equations cm be wztttm ir» differential afteraeor form 

&s 

(p^ - kg) s - < 2 p 4 - kg) pk| - akg - bkg 

g 

( 2 p - k 3 ) as 4 * <p - kg) y » k^ - s & 3 - bkg 



The characteristic equation o& this &yazam is 



<P^’ ** kg) <P^ ** k 5 ) 4- (2p 4 k 3 ) (2p - K^) « 0 
or 

4 . (4 «* kg - kg) 4 (kgkg • k^) « 0 
and the roots of the cheraeteristic equation ore give®, by 




4-kg-Ug 

^MlCMCtKUMID 
2 




(k 2 k 3 -kg) 



SO 



;t 11 

A stable solution is obtained if 



(4-Up-k ,f 2 

£-2- - k 2 i< 5 ik,70 



Substituting for the k'e from equations 8.7 through 8.11, 
the stability criterion is 

'fj, 2 - 



<4 + r/-*<4+A>- 



r; 



r: 



r; 



3 — 5- 

R T"i 



7 ° 

(8.14) 



Equations 8.12 and 8.13 can be solved to obtain a solu- 
tion in the form of constants plus pure sinusoids if the 
criterion of Eq« 8.14 is met. If the criterion is not met, 
the solution contains exponentially increasing sinusoids. 

Whether or not any significance can be attached to this 
criterion is not immediately apparent. If the solution resulting 
from a linearization of the equations about a point displaced 
from the Trojan point is to be considered better than the solu- 
tion of Chapter 5, it can be used only over © small segment of 
the orbit, and then a linearization in the vicinity of a net*? 
point must be accomplished. We cannot then categorically 
state that the stability of a specific periodic solution is 
significant if we consider only a small fraction of one period 
of the solution. 

On the other hand, it is felt that, if the vehicle passes 
through a region where the solution of its linearized equations 
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is unstable., there should be a tendency toward stability of 
the complete orbit. It seems that each time the vehicle 
passes through this region,, there would be & cusaslafcive ten- 
dency for it to be pushed further away from the Trojan point. 

Time was not available for a survey of the region in the 
vicinity of the Trojan point to determine where the criterion 
of Eq. 8.14 is set. Roughly it seems that it td.ll not be met 
for either x or y greater than about o.l. Xn any case, we 
should first have a better justification of the significance 
of such a survey. Although a time limit cm the present study 
precluded further investigation along these lines, it le felt 
that there may be something useful in this area. 
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C&t&FTElL 9 




At this point# the ignorance of the sutfeor o£ the infer i- 
cycles of classical celestial taoohsn ics is probably already 
apparel* 'Eh® literature has bam searched, but not thoroughly 
digested,.. If the method by which cur problem will fe© solved 
Kcsains hidden in this field, it can ssogft probably b© found in 
the work of Ausrctl Minfcnar (7) . 



,, The reeeafc Russian work on the subject has not hmn tho- 
roughly esantfxiQd by the author, partly because of tkm noa-svalXf 
feilifcy off fermslsfctas* The work of Mermen (S) is briefly di®~ 



cussed by X-eimnis (9) . This method seems to be &1 bq described 
by m Abstract off a paper by Schigolenr (10); 



”Tho action of ®a asteroid its calculated according to 
approximation formulas ..... -The *ppr©&&aaasi©» values desired 
are obtained by averting. Th® average method of the author# 
Ssowevar# differs ffsrcsa the vmll-k&om msfehod off ceLsselal 
aaechaiaice inasmuch as f;ha author does not average the poster 
function# but certain skillfully chosen parto of the differen- 
tial equations and thus obtains solvable differential equafeioag” 



This method was 



*5 Wi^Unrv /■?•• V* 






periodic solutions were indeed obtained# but not far enough to 
jt*dge the validity of these solutions . 
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iMinm io 

**'»**«;*»*, n >00. -wiJ'.sit' •-pvtjx 



co&c&ssxcfts 



7 



A8 stated isi Chapter X* wa do not have specific answers 
t:o the problem uzzdei consideration- if wa are asdcsd ,J Hos? 
c®sa the orbit©.! stability of a proposed Tzojm point satellite 
ba predicted?’ 5 , we can suggest the following approach; 

i.« Evaluate the expected guidance inaccuracy in terras 
of velocity md diapXaecsM&i error at the ‘Szojm poi**t„ Select 
Chs aiding point os suggested in Section 5*« to ®d»imi»e the 
else of the resulting orbit about the fmjsa point., based on 
tbs solution of the linearised ctpmfioai s* 

2* Compute the orbit around the Trojac point, using as 
incremental sreferenna orbit®, solutions of the equations of 
Chapter 8* Th® length of the tajejMeats used in the compu- 
tation will be picked te mini miss the cojsbiaatioa of truncation 
error and rounding error. Also m atiaapt should bo gaadn to 
predict ho# the sis® of this orror will increase with time* 

If the eenpatad orbit divarg©* from the Trojan point before fcrdr- 
co&put&tica error becomes bigniflcoi.it, tto question ta aaswer*td 

But a correct negative aucrer is much easier to obtain SSmk* 
a correct poative mmmT» if t£se emitted orbit diverges to 
the Trojan point only slightly, it will rot hr. «X#a» whether 
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this is du@ to m aefcuoX SasfcpMXifcy of the arbifc, or a result 
of the accumulated eos^sutati ion error* 

3* Fust tlier work on the orbital stability of Trojan 
^ > satellites should probably be directed toward the problem of 
.dentlfyiag and reducing the effect of computational errors* 



it 



Also it is felt that further investigation should be raade 
of the significance of the region in which stable solutions 
of the linearised equations are obtained* as discussed in 
Chapter 8* 
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The following calculation;' w©r e aada to permit calcu- 
lation of a collation o.f the linearised equations dirastfly 
from specified initial eondifcioas* For ni*aar.‘cal worh, a 
lvalue of s of exactly 0 , 0,12 13 was need. AX&iirugh, the 
physical quantity in the earth '&oon system i. r only knarni 
to within about 0.3 par sent* the calouto'&ach wra carried 
to ton significant figures to permit uw of the results ia 
subsequent nuaarieal integre. -low, ccsaput/ticuo, if desired. 
Tte# result* are then GOficlnfeent with a linearised model 
of the system uifch tho mma ratio upm: ttlmd by € m 0.01213 
exactly. 



A 



& i 



0.98787 



H 1-272 ( / - e ) « 0 , 3214735232 



pf ••*• % <1 h 0.8224735232) => 0,91X2367616 



p| » h (1-0.8224735232) 



0 *0087632386 



fS3 



p x « 0.9545872205 



p 2 = 0.2979316002 



«5T» 

o^lo 

«“ (1-26) = 1,267523442 
4 



Substituting initial conditions s < 
y (0) into aquations 5.1 sad 5.2 


o 

V 

O 

*»• 

/% 

o 


36 (0) » B-£ 4* 


<A,1) 


y (0) « B « 4* £>,, 


(A. 2) 


K (0) 33 4* PpG-; 


(A. 3} 


7 (0) « PjAg + P 2 C 2 


(A.4) 



Combining equations 5.3 ®d 5,4 we can get' A^ in 
tamo of B-jj and Bg 
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«* 



a /% 9 
2p x A x « ~£LjU (l»2e) B x 4* <pj 4- 

A . 


8, . 

V J, 2 


(A.5) 


Sis&teXy fe/a ©amfitea® 5.5 smd 5*6 






*P 2 % CT <M e ) % * <pf * j) ® z 


(A. 6) 


Fma &/&&&&&& A.1 a; 4 A*- 2 






o js (0’i - Bj 




(A. 7) 


». t » y 0) - 3g, 




(AeS) 



9ab«fc£tut;ins ©qmfcisma A* 7 &&c1 A, 8 isfco Mq> A* 6 

a: (2-^S) jx (0) - B^j 4* (pg «> '«) |y(0) «&g | 

(A.9J 

*Sttlfciplyfcfj 8q. /.*3 by 2 m3 su&ee&Cutiqg 3teee ©quaftiona A#5 
A*9 

r=** ' - 

2s <o) ■•> -2^ js^ + ss(o)«>B^j-t' Cr|4)^Kpl 4* j) 

{A# 10) 



2 zz (G) & s?£ 0 ) 

4 






•Kp? •* ?|) %Kp|* *^)y(o) 

4 

<A.ll) 



^ S?*£0^' 

6 



2 2 
H p 2 




y(0) V -4— n 3 (0) 

pi ~ p a 

<A,12) 



Subasl&afci&g smst sgtcol values 

B 2 »1 ♦541111544 &<0) -2»349572&77 y(0) * ?.431<&SS219 k (0) 

* 

(A* 13) 



GossM&issg equations 3.3 a&d 5-4 to el£c&&a£@ 
of £±2* Gim get Ag in fcesass of B-, sad Bg 

O •». O \Aa 

52 •* (p« 4 1 ***) Bfl rj -'■■**■£> (l-2«) B>j 

«fi» «£• 4»* ^ *«• 



instead 



(A»14) 



Similarly £soja ^u&ti&s& 5*3 ©ad 5*6 



2p 2 C 2 » » <p§ V ~! »3_ 



Cl *wxn.'3irt 



<Ms) |> 2 (A. 15) 
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Substituting equations A./ and A.ii into Ecj. A. 15 
2© 



P2«2 * " < 4 + &-> '• “J- « ' *J 



(A. 16) 

Multiplying Eq. A. A by 2 and substituting from equations 
A. 14 and A. 16 

2y (0) - - (pf + B x • (3.-26) Jb 2 + y(0) - B 2 J 



-(p| 4- j^°) ” B lJ 



(A. 17) 



2y (v>) ® - 



i — ‘ 

(1-2C) y(o) - (pf"p|) B x - (p| + j) s(0) 



(A. 13) 



B, 



2 3 
Pp + 7 

— P — V~ x (°) - 
p - <iZ 
1 1 * 2 



3^ 

zzza=2£i y(U ),_L, 
2 2 yv _2 2 
P 1 " p 2 



p l“ p 2 



(A. 19) 



y<0) 



Substituting nuraerical values 

o -1.019305762 k( 0) -1.541111544 y(u) •• 2.431689219 y(u) 



(A. 20) 



can nmi be found from Eq» A* 5 

2 9 






A^ 83 



i »?<* ) 



., t Pi + 4 

1 2p, 



B. 



'2 



(A. 21) 



A x o 0.6639118009 B x -f 1,653813473 B 2 (A. 22) 

Substituting from equations A. 20 and A. 13 

A x “ 3*228854338 x(0) - 5*731588091 y(0) 

•3-4.026423771 i(0) - 1.614427169 y(0) 

(A. 23) 



Similarly from Eq. A. 14 



A 2 ” 



2 . 3 
Pi 4 * 4 

“ ~*^1 



B 1 ~ 



3 fl 

~IZJ2z£l B 



2Pi 



2 



(A. 24) 



A 2 *= -0.8701335645 B x - 0.6639118009 B £ 



(A. 25) 



Substituting from equations A. 20 end A. 13 
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Ag » 1.910529363 n<0) * 3.228S54338 y(0) 

-1.614427169 x (0) + 2.11589440S y (0) (A. 26) 

Substituting Eq. A, 20 into Eq. 4.7 

• 2.019805762 x<0) 4- 1.541111544 y(0) * 2.431689219 y (0) 

(A* 27) 



Substituting Eq. A«13 into 2q* A>8 

D ? « 1.541111544 n (0) -f 3.843572677 y (0) - 2.431689219 k <0) 

(A. 28) 




can be fennel fxeta Eq. A. 6 



C 



1 




( 1 - 2 ®) 




D 3 * 




9 

4 





(A. 29) 



C 1 *= 2.127205440 * 3.923000296 B fJ , 



(A. 30) 
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. 














Substituting from equations A. 27 and A. 2d 



C x ® 10.345405o7 x<0) -5- Id. 364234/6 y(0) 

“9 .544380904 ::(o) + 5.172/02535 y(o) 

(A. 31) 



Similarly from Eq. A, 15 



^2 a ~ 



2 , 3 
p 2* 4 

'~wr~ 



D 



1 



3 3 

-Z3^LO=2£l D , 

2p 2 






C 2 " 



- 1. 4o7643965 D-, - 2.127205440 D 



2 



(A. 32) 
(A. 33) 



Substituting from equations A. 27 and A. 28 



02 ”=“ 6.121428251 2 c(o) - lo.345405o7 y(o) 

* 5.1/2702535 x (o) - 3.422952654 y(o) 

(A. 34) 
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Givem specified 'Initial euanditiens* the constants 
of integration eon now fee calculated from equations A* X3 5 
A*2Qj, A, 25 3 A*27* A* 28. A«JX and A*34* The so equa- 

tions war© crosschecked roxaerieallj by other fonoolas and 
found to bo consistent* 



The substitution of aml« into this aquation tells 
us little about mx&i: ©ost of: m. orbit we would gst with 
slightly different Initial /.©adifcions* as does the saetbod 
of Section 5*5. Bouever* :/or calenlarion of a specific 
©Kbit, the a^serica,1. w/dt ‘Is aufth simpler* If this solo," 
fcioa war© required an a ruferenoe orMt for a digital 
coraputafeion of fchs orbit* the abo'*?a €,<jnatioEs would fee used. 
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APPEOSXK B 



imwiicAh camwtm of tm 
— — — — - — >«««.■ • — 



MomuM mm^im 



To defce; raise scrs® m^sricaX values of u a&d u% 
orations 4.5 easd 4.20 xs«*e feranafoasrasad to polsr coosrdi- 
i'iiter. by tfre equations 




y « jp sin 



The equations &h@n become 



u *> 



e ,o cos % ( j? *cln^ 6 3 f fife ^ * eoa^$ hh jj cos P •*■!} 

4y> sfe®* J) 2 csos0*ycas£+ 1} 



***■5 



♦ e ( ® *«iA 6p 8ta©+J>® ooe*d-j> «og^ 1)” 

<S.l) 

^ J7~ ^ ^ 

u* « ~ P^coe 2 #* <X*2«)P 2 8in & sos$ * - p 4 sj«? £ 
2 S •’ 4 ' $ * 

(B«2) 
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V 


























Using she identities 




4* 




1 




the aquations 



£02? n sand u* hQC-Gfi.&s.i 



u 







-60°) 4- X I 



■4/{ 






pcm 



( 0 -60°) 





<B*3) 



$ 

» *3 









{B»4) 



Using a value fos? ® of 0 ,0X2X3., file cguG&St&ie® u ©ad u 5 
war© cw&lcuXs&ed for 0*1 isicaiOEKrafcfi of te 0 f:o Q*S mid feoc 



30 ° iA&sm&a&Q qS & * 
Tho ae&idu al 



%h® results are Xictcd ie 5P«ibl 
u" «* « - a* at fche#e assists 1® 



.© B*l 
shawa 



in Fig* 
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Distances from the trojan i^issc in 0.1 incraes&nfcs 



^ A 

Indicated values are u r: IQr 



Fig. B. 1 % Appr ©xissst© Values of the Residual Jacobian Function 



68 



TABLE B .l 

TABUL ATION O? JACOBIAH FUNSTIOHS 
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